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Light deflection in the post-linear gravitational field of two bounded point-like masses is treated. 
Both the light source and the observer are assumed to be located at infinity in an asymptotically 
flat space. The equations of light propagation are explicitly integrated to the second order in G/c 2 . 
Some of the integrals are evaluated by making use of an expansion in powers of the ratio of the 
relative separation distance to the impact parameter (ri2/£). A discussion of which orders must be 
retained to be consistent with the expansion in terms of G/c 2 is given. It is shown that the expression 
obtained in this paper for the angle of light deflection is fully equivalent to the expression obtained 
by Kopeikin and Schafer up to the order given there. The deflection angle takes a particularly simple 
form for a light ray originally propagating orthogonal to the orbital plane of a binary with equal 

in ■ 

masses. Application of the formulae for the deflection angle to the double pulsar PSR J0737-3039 
for an impact parameter five times greater than the relative separation distance of the binary's 
^— ^ ' components shows that the corrections to the Epstein-Shapiro light deflection angle of about 10 -6 

arcsec lie between 10 -7 and 10 -8 arcsec. 

>v 

j3 . PACS numbers: 04.20.Cv, 04.25.-g, 04.25.Nx 
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■ I. INTRODUCTION 

<n : 

Light deflection by a gravitational field pj is one of the observational cornerstones of general relativity. The 
observational confirmation of Einstein's prediction that light would be deflected by the gravitational field of the Sun 
brought general relativity to the attention of the general public in the 1920's. Today technology has reached a 
level at which the extremely high precision of current ground-based radio interferometric astronomical observations 
approaches 1 /xarcsec and within the next decade the accuracy of space-based astrometric positional observations is 
also expected to reach this accuracy. At this level of accuracy we can no longer treat the gravitational field of a system 
of moving bodies as static and spherically symmetric. This fact is one of the principal reasons for the necessity of 
a more accurate solution to the problem of the propagation of electromagnetic waves in non-stationary gravitational 
fields of celestial bodies. To reach the accuracy of 1 /iarcsec, many subtle relativistic effects must be taken into account 
in the treatment of light propagation in non-stationary gravitational fields of moving bodies. One of the most intricate 
| problems is the computation of the effects of translational motion of the gravitatingbodies on light propagation. 

This question was treated for the first time by Hellings in 1986 0- In 1989 KlionerQ solved the problem completely 
to the first post-Newtonian (1PN) order (i.e. to the order 1/c 2 ) for the case of bodies moving with constant velocity. 
The complete solution of the problem for arbitrarily moving bodies in the first post-Minkowskian approximation 
(linear in the gravitational constant G) was found by Kopeikin and Schafer in 1999 They succeeded in integrating 
analytically the post-Minkowskian equations of light propagation in the field of arbitrarily moving masses. In Rcf. 

0, Le Poncin-Lafitte et al. have recently developed an alternative approach to the problem of light deflection and 
time/frequency transfer in post-Minkowskian gravitational fields based on an expansion of the Synge world function 
for null geodesies. In that paper the world function and time transfer function were computed for a static spherically 
symmetric body to the second post-Minkowskian approximation. 

In this paper we treat light deflection in the post-linear gravitational field of two bounded point-like masses (binary 
system). We assume that the light source as well the observer are located at infinity in an asymptotically flat space. 
To compute the light deflection we integrate the equations of light propagation explicitly to the second order in G/c 2 , 

1. e. to the order G 2 /c 4 . 

The assumption that the gravitational field is weak along the light path allows us to consider the metric as a 
perturbation of a flat metric represented by a power series in the gravitational constant G 

oo 

g MV [^, G] = 4°J + ]T G»g$ {x°), (1) 
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For the same reason we can consider the light trajectory as a perturbation of its trajectory in flat space (a straight 
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line) represented by a power series in G 

oo 

^) = ? ( o)W+E G "4)W- ( 2 ) 

71=1 

It follows from Eq. @ that the vector tangent to the light trajectory takes the form 

f (t) = ^ = l m + Y^G n 5l {n) {t), (3) 

n—l 

where £(q) is the constant vector tangent to the unperturbed light trajectory. 

In order to obtain the post-linear equations for light propagation we introduce into the differential equations for the 
null geodesies the metric as given in Q and the expression for the tangent vector l(t) as given in ©. As a result we 
get a set of ordinary coupled differential equations of first order for the perturbation terms 5l^(t). Each term <%(„) is 

given in the form of a line integral along a straight line in the fictitious metric gffi , i.e. along the original unperturbed 
light trajectory. We get the post-linear light deflection to the order G 2 /c 4 after computing the perturbation terms 

(5l(i)(i), 5l(2)(t) and the corrections arising from introducing the linear perturbation of the light trajectory 5z^(t), 
the motion of the masses and the shift of the lPN-centre of mass with respect to the Newtonian centre of mass in the 
expression for linear light deflection. The final result we obtain is the expression for light deflection in the post-linear 
gravitational field of two bounded point-like masses. The deflection angle takes a particularly simple form for a light 
ray originally propagating orthogonal to the orbital plane of a binary with equal masses. 

This paper is organized as follows. In Section [H] we derive the post-linear light propagation equations. We set 
out an approximation scheme to integrate these equations. The deflection angle as a function of the perturbations 
of the vector tangent to the light ray is introduced. In Section IlIII the post-linear and linear metric for two bounded 
point-like masses in harmonic coordinates are given. The coordinate frame is chosen so that the lPN-centre of mass 
is at rest at the origin. In Section Hvl we compute the perturbation of the vector tangent to the unperturbed light ray 
and the corresponding light deflection in the linear gravitational field. In Section we compute the light deflection 
in the post-linear gravitational field. To facilitate the computations we separate the light deflection terms which are 
functions of the post-linear metric coefficients from the terms that are functions of the linear metric coefficients and 
the perturbations of the first order in G of the vector tangent to the unperturbed light ray. The resulting integrals are 
given in Appendices IA1 and IB1 In Section fVTl we calculate the additional linear and post-linear light deflection terms 
arising from the introduction of the motion of the masses into the expression for the linear perturbation. In Section 
I VIII we compute the corrections to the post-linear light deflection arising from the introduction of the linear perturbed 
light trajectory into the expression for the linear light deflection. The resulting integrals are given in Appendix IU1 
Section lVllIl is devoted to the computation of the corrections to the linear and post-linear light deflection arising from 
the introduction of the shift of the lPN-centre of mass with respect to the Newtonian centre of mass into the expression 
for the linear light deflection. The resulting expressions for the total linear perturbation, linear and post-linear light 
deflection by the gravitational field of two bounded point-like masses are given in an explicit form in Sees IIXI and Q 
In Section IXll we present our results and give in an explicit form the light deflection expression for some simple cases. 
The derived formulae for the angle of light deflection are applied to the double pulsar PSR J0737-3039. In Appendix 
[0]we compute the linear light deflection terms arising from the acceleration terms in the metric coefficients h^J and 
h^q . In Appendix |E] we show that the expression for the linear light deflection given in this paper is fully equivalent 
to the expression for the light deflection computed by Kopeikin and Schafer |5j in the event that the velocities of the 
masses are small with respect to the velocity of light and the retarded times in the expression given by Kopeikin and 
Schafer are close to the time of closest approach of the unperturbed light ray to the origin of the coordinate system. 
Finally, Section IXIIl is devoted to the discussion of the results. 

Notation 

Let us summarize the notation and symbols used in this paper: 

1. G is the Newtonian constant of gravitation; 

2. c is the velocity of light; 

3. the Greek indices a, f3, 7..., are space-time indices and run from to 3; 
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4. the Latin indices i, j, k,..., are spatial indices and run from 1 to 3; 

5. g^v is a metric tensor of curved, four-dimensional space-time, depending on spatial coordinates and time; 

6. the signature adopted for g M „ is ( — h ++); 

7. we suppose that space-time is covered by a harmonic coordinate system system (a; M ) = (x , x 1 ), where x° = ct, 
t being the time coordinate; 

8. the three-dimensional quantities (3- vectors) are denoted by a = a 1 ; 

9. the three-dimensional unit vector in the direction of a is denoted by e a = e* ; 

10. the Latin indices are lowered and raised by means of the unit matrix Sij = 5 lJ = diag(l, 1, 1); 

11. by ,o- we denote the partial derivative with respect to the coordinate x a ; 

12. the scalar product of any two 3-vectors a and b with respect to the Euclidean metric 5{j is denoted by a • b and 
can be computed as a ■ b = SijtfV = a l 6 8 ; 

r iV2 

13. the Euclidean norm of a 3-vector a is denoted by a = \a\ and can be computed as a = o mn a m a r ' 

14. by Z(q) we denote the vector tangent to the unperturbed light ray z(t) and the unit vector e*( ) is defined by 

6(0) = £(0)/R(0)|- 

II. LIGHT PROPAGATION AND LIGHT DEFLECTION IN THE POST-LINEAR GRAVITATIONAL 

FIELD 

A. The light propagation equation 

In this paper we calculate the light deflection in the post-linear gravitational held of two bounded masses for the case 
when the impact parameter |£ | is much larger (5 times or more) than the distance r\i between the two accelerating 
masses, so that we can suppose that the gravitational held is weak along the light path. 

For weak gravitational holds we can assume that the light propagation is very well governed by the laws of geometric 
optics, whereby light rays (photons) move in curved space-time along null geodesies. The equations of null geodesies 
with the time coordinate as a parameter are given by (e.g. see 0) 

—+Yy^ = c -^i a m\ (4) 

where 

r p<r = ^5 mA [<?pA, ct + g*\, P - g P a.x] (5) 

are the Christoffel symbols of the second kind and = dz^/dt denotes the 4-vector P = (c,l l ). Here, it is important 
to recall that — (z°,z l ), where z° — ct. Notice that l M is not exactly a 4-vector because we differentiate with 
respect to the time coordinate t. So I 1 * is a 4-vector up to a factor. The space component of Z M given by P = dz l /dt 
is the 3-vector tangent to the light ray z l (t). In the present case of null geodesies, Z M has to fulfil the condition 

l 2 =g^[z°,z\t),G}l^r = 0. (6) 

Now we consider a light ray z l (t) that is propagating in a curved space-time g^z , z l (t) 7 G] with the signature 
( — h ++). If the gravitational held is weak, we can write the fundamental metric tensor g^z , z l (t), G] as a power 
series in the gravitational constant G 

oo 

9iiV \z\ z\t),G] =V^ + Y1 ^"V, zi (t), G], (7) 

71=1 
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where rj^ v is the Minkowski metric and hffi [z°, z l (t), G] is a perturbation of the order n in the gravitational constant G 
equivalent to G n g'${z , z l (t)) of Eq. (JIJ (physically, this means an expansion in the dimensionless parameter Gm/c 2 d 
which is usually very small, d being the characteristic length of the problem and m a characteristic mass). 

In order to obtain from equation (J2J the equations of light propagation for the metric given in 10 we substitute 

the Christoffel symbols into the equation To save writing we denote the metric coefficients h\}q [z°,z*(t),G], 
h p 2 q [z°, z l (t), G] by h^q and h p 2 q . Then the resulting equation of light propagation to the second order in G/c 2 is 
given by 

dl% 1 „2i,(l) 2,(1) ,(1) ,m , mi ,(1) ,m h W 



dt 2 



c n 00,i c n 0i,0 L/i 0i,m l + L "Om,i' L ' t mi,0 ( n mi,n l 1 



,\W /mm _ 1-1.(1) ,i _ ,(1) .fc.j /l X ,(l) -1,(1) \/m,p/i 

^ 2 mn,i l 1 2 00, l "-OO./c 1 1 ^2 mp,Q L n op,mJ L 1 1 

+ - c 2 h {2) - - c 2 h^ ik h {1) - h (2) i k r - ( h {2) --h {2) )rr 

' 2 00, i 2 00, fc ""00, fc 1 y L mi,n ^ rnn,ij l 1 

+ h {1)ik (h% n - \h { l, k )i™i n - h®h<»,M, (8) 

where by q and , we denote d/dz° and djdz % respectively. To calculate the light deflection we need to solve equation 
© for Z*. In order to solve this complicated, non- linear differential equation, we turn to approximation techniques. 

B. The Approximation Scheme 

We can write the 3- vector l l (t) as 

= *{o) + !>(»)(*)- ( 9 ) 

n=l 

where ZLs denotes the constant incoming tangent vector oo) and Sl^(t) the perturbation of the constant tangent 
vector Z^ of order n in G equivalent to G n 8l^(t). 

After introducing the expression for l l (t) given by Q into the equation © we obtain differential equations for the 
perturbations Sl 1 ^ and 8l^ 2 y These are given by: 

<^(l) — 1 r 2 hW -r 2 h W -rhW l m ^rh {1) l m -r/i (1) l m - /i (1) l m /" 
— 2 00,4 6 "Oi.O c ' t 0i,m t (0) + L,l 0m,i l (0) c ' t mi,0 ( (0) n mi,n l {0) '(0) 

, lft(l) /m in _1„J,(1) /» _/,(!) /ft 7' _|_ fl^-l 7, W -r~l/) (1) \im IP n (-\C\\ 

+ 2 m ™. i (°) (°) 2 00 '° (o) n, oo,fc t (o)'(o) + ^ ™p,o c n 0p,mj l (o) l (oy(o) \ w ) 

and 

- 1 r 2i ) (2) _ 1 2, (i)ifc, (l) _ h m jk ,i _( h m _V 2) "l/ m 7™ 

^ — 2 00,i 2 00, fc /t 00,fc ,, (0)''(0) ^"mt,n 2 mn ,i ) (°) (°) 

1 h (l)ik( h (V 1 l(1) \/mm 1,(1)7,(1) /fc /i 

+ " ^ n mfe,n _ 2 m "' fc / (°) (°) 00 n 00,fc t (0) < (0) 

„/,(!) s/m 1 „/,(!) X/m ,.7,(1) A7' 
_ c "0i,m°'(l) + cn 0m,i° l (l) Cll mi,0° l (l) 

1) JJ/1 

1 



1,(1) Jf/m m i(l) im r,n , , (1) r,m in 
"mi.n"^!)'^) ~~ "iih,h 1 (0) <,1 (1) + n mn,i 0L (1) L (0) 



lr>7)( 1 ) A7* Ti^ 1 ) A7 fc 7* 7 ) ( 1 ) 7 fe A7* 
2 C "oo,o 0( (i) ~ "oo,fc ot (i) ( (o) ~ "oo,fc'(o) 0f (i) 



+ c 1 ' l Lp,0^(l)^fo)^(0) _ c 1 hap,mM™) ^(0)^(0) _ c I ' l ip,m'(0)^fl)'(0) 

+ (2° ^"ip.O" I ' l 0p!m) '(O)^fo)^(l)- (^) 
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In order to calculate the perturbations 51^ (t) and 8l 1 ^ (t) we have to integrate equations (|10|) and along the 
light ray trajectory to the appropriate order. 

Before performing the integration it is convenient to introduce a new independent parameter r along the photon's 
trajectory as defined by Kopeikin and Schafer ||. The relationship between the parameter t and the time coordinate 
t is given by 



t-t* 



(12) 



where t* is the time of closest approach of the unperturbed trajectory of the photon to the origin in an asymptotically 
flat harmonic coordinate system. Then the equation of the unperturbed light ray can be represented as 



z i (r)=Tl\ 0) +S 



(13) 



where is a vector directed from the origin of the coordinate system towards the point of closest approach. The 
vector £ l is often called the impact parameter and is orthogonal to the vector l^ y The distance r(r) = |z(r)|, of the 
photon from the origin of the coordinate system reads 



r(r) = y/c 2 r 2 + £ 2 



(14) 



It follows from equation l|12|) that the differential identity dt = dr is valid, so that we can always replace the 
integration along the unperturbed light ray with respect to t by the integration with respect to the variable r. 
Then the resulting expression for 81)^ is given by 



ch {1) 

c,l 0i 



da K 



(1) n I 

oo,o l (o) l(-0 "t" 



da I 



(0)'(0) 



lm 

n mi L (0) 



"00 '(0) 



mp,0 



L n 0p,m 



l (Q) l(-) 



(15) 



On the right-hand side of equation 115|) after evaluating the partial derivatives of the metric coefficients with respect 
to the photon's coordinates (i.e. (z°,z l (t))), we replace in the integrals the photon trajectory by its unperturbed 
approximation z l (cr) unport . = al 1 ^ + £ l and the time coordinate z° by a + t* . In this paper we denote this operation 
by the symbol |(— >)• Then we perform the integration with respect to a. 

After substituting the expression obtained for 8l 1 ^ into equation l|llfl we can integrate it to get 8l l , 2 y To calculate 
the perturbation 81^ we separate the part of 8l 1 ^ that depends on the post-linear metric coefficients from the part 
that depends on the linear metric coefficients. We denote these parts of Sl 1 ^ by <W 2 \j and <W 2 )n respectively. As in 
the case of equation l|15H we replace the photon trajectory by its unperturbed approximation and the time coordinate 
z° by a + 1* after evaluating the partial derivatives of the metric coefficients with respect to the photon coordinates. 
The expressions for Sl 1 ^ and <W 2 N n are given by 



<4)i(t) = 



da 



rl 



2 C "•oo,i "oo.fc'OT'Co) 



(2) 



da 



im in I 
'(O)'(O)K^) 



(16) 



and 



<«(2)IlM 



da 



3t(l)ifc fc (l) , 1,(1)1,(1) ,k ,i 
iii "•nn h i "fin "nn fc'fnWn 



da 

da 
da 

i 

da 



1 

2 



^(Cn-^^UJ^)%lH) 



^(i) 

"0m,i "0i,m "'mi,0 



(i) 



k-) 



i 



-ch^Sl^ia) + h$ ik 5lfa(*)r m + K^Sl^ia) | H) 



,(1) JTlk 



(1) rfe A ;i 
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+ c 



da 



da 



Z (o) ? (o)^(i)( CT )l(-)- 



(17) 



C. Light deflection 

The dimensionless vector a 1 ,^ of order n in G describing the angle of total deflection of the light ray measured at 
the point of observation and calculated with respect to the vector l 1 ^ (see 0) is given by 

a\ n) (T)=P>^f- : (18) 



where 51^ is the perturbation of the constant tangent vector of order n in G. Here, P l q — S q — e^e( ) g is the 
projection tensor onto the plane orthogonal to the vector K Q y In the case of light rays (photons) | Z > I = c - 



III. THE POST-LINEAR GRAVITATIONAL FIELD OF TWO BOUNDED POINT-LIKE MASSES 

In the computation of the metric generated by a system of two bounded point-like masses we distinguish between 
3 zones: the near zone, the intermediate zone and the wave zone or far zone. 

In Refs. and it was shown that leading order terms for the effect of light deflection in the case of small impact 
parameter depend neither on the radiative part (~ l/£) of the gravitational field nor on the intermediate (~ l/£ 2 ) 
zone terms. The main effect rather comes from the near-zone (~ l/£ 3 ) terms. Taking into account this property 
of strong suppression of the influence of gravitational waves on the light propagation, we can assume in the present 
work that the light deflection in the post-linear gravitational field of two point-like masses is mainly determined by 
the near-zone metric. 



A. The metric in the near zone 



In Ref. 10], Blanchet et al. calculated the conservative 2 PN harmonic-coordinate metric for the near-zone of a 
system of two bounded point-like masses as a function of the position z and of the positions and velocities of the 
masses x a (t) and v a (t) respectively, with a = 1, 2. We shall use their metric in this paper. The post-linear metric for 
two bounded point-like masses, at the 2PN-order, is given by 



u(2) 

'on 



G 2 m\ 



G m,\m 2 



r%T2 



ri 

zr 12 



2r2r 3 



12 



2r 2 ri2 





6 pq 


rG 2 m 




- r\ 


G 2 m\ . 


i a 




'i 





-i + G 2 m 1 m 2 l 



r\r 2 



2r 2 r\ 2 



2riri2 r 12 S 



G2m l ^P r q AC 2 m, mnn p n q ( 1 4- 1 \ , 4G 2 miTO 2 ( (p a) ( P a)' 
— n 1 n 1 - 4Ct mim 2 n 12 n 12 ^— + — — — J H — yn 1 n 2 + Zn 1 n 12 



(19) 



(20) 



where the symbol (1 «-> 2) refers to the preceding term in braces but with the labels 1 and 2 exchanged; by S we 



denote S = ri 



r 2 



ri2, where n 



XI 



(*)l> r 2 



n p 2 are unit vectors defined by n p = r p /r\, n 2 = r 2 /r 2 and n p 2 = r p 2 /r\ 2 . 



x 2 {t)\ and ri2 = \x\{t) — x 2 (t)\. The vectors n p , n 2 and 
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In our computations we also need a part of the linear gravitational field of two bounded point-like masses. The 
part that is relevant to our calculation is given by 



x (1) - 



2 2 

„ G \ ^ m a G v— * rn a r 

a— I a—1 



2 1 4^ 2 



m 



(21) 
(22) 
(23) 



a—1 a—1 

where denotes the velocity of the mass m a and is the unit vector defined by n v a — r v a jr a . 

Here, it is worthwhile to point out that the parts of the linear gravitational field in h^j and h^q that contain the 
accelerations of the masses were introduced into the part of the gravitational field quadratic in G after substituting 
the accelerations by explicit functions of the coordinate positions of the masses by means of the Newtonian equations 
of motion. 



B. The barycentric coordinate system 



We use a harmonic coordinate system, the origin of which coincides with the lPN-centre of mass. Using the 1PN- 
accurate centre of mass theorem of Ref. jllf , we can express the individual centre of mass frame positions of the two 
masses in terms of the relative position ?xi = x\ — x,2 and the relative velocity v\2 = v\ — V2 
as 



where X\, X2 and c\pm are given by 



Here, we have introduced 



Xi 
X2 



X 2 + £ipn 



ri2, 



X\ + -rClPN 



1 M 



Xn 



£ipn 



m 2 
v{m\ 



m 2 



2M 



GM 



12 



ri2 



(24) 
(25) 

(26) 
(27) 

(28) 



M = mi + m 2 , V12 = Ivu] 



(29) 



and 



TO1TO2 
A/ 2 ' 



(30) 



Here, it is important to remark that in our computation of the post-linear light deflection up to the order G 2 /c 4 , we 



need only to consider the lPN-corrections to the Newtonian centre of mass, because, as we shall see in Sec lVIHl the 
2PN-corrections to the Newtonian centre of mass are related to post-linear light deflection terms of order higher than 
G 2 /c 4 . 



IV. LIGHT DEFLECTION IN THE LINEAR GRAVITATIONAL FIELD OF TWO BOUNDED 

POINT-LIKE MASSES 

In this section we compute the perturbation term 5l' 1 ^ (r) and the corresponding angle of light deflection for an 
observer situated at infinity in an asymptotically flat space. To compute the perturbation term S1 1 (1 Jt) we have 
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to introduce the linear metric coefficients given by Eqs (|21(l - (|23|l into Eq. (|15|l . After introducing the linear metric 
coefficients into Eq. 1)150. we find 



2 

dl\ 1) (r) = -2Gj2m a f da 1 [*« - x*(t)] | H) 

G , 2 . , /" T 1 
+ 4 -^m Q / <&r— (e (0) ■#«(<))[>' -a£(*)] Ih) 

3 k(t) + (e (0) • i(t)) 2 ] [z 4 - xl(t)] | M 

a 

3 §E m ° / ^i(^-^(*)) 2 [^-4(*)]l(-) 



2 §E m ° 

° 0=1 ' v 

2 



- 4 §E m ° 

° 0=1 ' ' 

C 0=1 



dcr 13^(0) ■«<»(*)) ccr- (e (0) -f a (t)) /( 0) | ( ^) 



(e( ) ■Va{t))ca+ £ ■ Vait) - X a (t) ■ V a (t) ><(i)|(^) 



a— 1 a— 1 a— I 



2 §E m °- 



-2f«2(*)ho)> 



(31) 



where r a = z - f a (t) and r a = \f a \. 

Because the linear metric coefficients are functions of the positions and velocities of the masses x a (t) and v a (t) 
respectively, the expression for 51 1 ^(t) given in Eq. 1)31(1 is a function of these quantities. This means that we have 
to take into account the motion of the masses when we are going to compute the integrals in Eq. H31[). Considering 
that the influence of the gravitational field on the light propagation is strongest near the barycentre of the binary and 
that the velocities of the masses are small with respect to the velocity of light, we are allowed to make the following 
approximations: 

1. We may assume that the linear gravitational field is determined by the positions and velocities of the masses 
taken at the time of closest approach (t = t*) of the unperturbed light ray to the barycentre of the binary (i.e. 
to the origin of the asymptotically flat harmonic coordinate system). The expression, resulting from Eq. I|31|) 
after setting t = t* for the positions and velocities and computing the integrals, is denoted by 51 (^(t); 

2. We treat the effect of the motion of the masses on light propagation as a correction to the expression of SI^^t), 
which we denote by <5^ 1 ^ II (r). We shall compute this correction in Section IVll 

After fixing the values of the quantities x a (t) and v a (t) to x a (t*) and v a (t*) and evaluating the integrals in Eq. H31fl . 
we find 



^ 1)I (r) = -2^m a So[f -<(0]-§E TO «{ 2Ai + r}' 



(0) 



C C T a 

a— 1 a— 1 

G 2 ^ r 1 1 

^L m »( e 1»] ,i? «( i, ))(A + -f(o) 

a— 1 a 
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G 



|6 ? ^m o (e (0) -i(f)) i-v a {t*)-x a {t*)-v a {t*) F a3 [? -<(**)] 



G 



0=1 
2 



2 -^E m °^*) s 4e -<(**)] 



G 



a=l 
2 



+ 3 ^E m « f #«(**) -*»(**) •*„(**) F a4 [e - x \{t 



0=1 
2 



2^a + — 



+ 2 



_ -3 E TO <M -ir «(**)) 
0=1 I 
2 

+ 6 - m.(e"(o) • $.(**)) • tT„(t*) - 2„(t*) ■ 

a— 1 

+ ^"» o (e (0) .i(f)) 2 {3F ol - 2A a )l\ 0) 

a— 1 

- ? E ™.«2(n{2 a, + f }iu + 2 ^ E ^(r- unfi 

a— 1 a— 1 u 

2 

+ 3 -4 E m « [£• - • * r «(**)l 

a— 1 

where the functions A a , B ai F a i 7 F a 2, F a 3, and F a 4 are given by 

r2(0,r) + (e (0) -x a (t*))( C r + r a ) 
(e (0 ) -x a (t*)) +cr + r a 



e-tr„(t*)-f„(t*)-tr„(t*)Js„|«i(t*) 

^a2^( ) 



(0) 



1 








1 




r a R a 




1 


/ 



3«| 

+ 2(e (0) .x a (t*)) 



CT 



2r*(0,t*) - (g^o) ■:?„(**)) (3r2(0,f) - (e (0) • x a (t* )) 2 ) r a 
3r 4 (0,t*) - (e (0) • £ a (**))(3r 2 (0,i*) - (e (0) • x a (t*)) 2 ) r a 
3^(0,i*) + 3(e (0) -x a (t*)) 2 r2(0,i*) - (e (0) • f„(t* )) (3 r 2 (0, f) - (e (0) • x a (t*)) 2 ) r a 
3r 2 (0,t*)-(e (0) -f a (t*)) 2 



2 2 

C T 



+ (e (0 ) -x a {t*)) 
1 f 



c 3 r 3 



a2 



3r 3 i? 2 



2(e (0) -f a (OK(0,O + r 2 (0,r) + (e (0) -f Q (i*)) 2 r 2 (0,i*)r a 



- 2 (e (0) • £„(**)) (r 2 (0, f) + (e (0) • x a {t*)f) r a - 3 (e (0) • x a (t*)) r 2 (0, t*) 
+ 



CT 



+ 



r 2 (0, f) + (e (0) • x a (t* )) 2 r a - 3 (e (0) • £„(**)) 
r 2 a (0,t*) + (e {oy x a (t*)) 2 



2 2 

C T 



c 3 r 3 



F a3 = ^j-^<! - r 2 (0,r) [r 2 (0,O + (e (0) • x a (r)) 2 - 2 (e (0) • f (t*))r„ 



+ (e (0) -x a (t*)) 3r 2 (0,r)+3(e (0) •£«(**)) - 4(e (0) -x a (t*))r a 



CT 



(32) 

(33) 
(34) 



(35) 



(36) 



10 



+ 2 (e (0) • x a (t*))[r a -3(e (0) • 2 a (t*))\ c 2 t 2 + 2 (e (0) • x a (t*))c 3 r 3 j , 
Fa4 = R~RM (e " (0) ' fa(r)) I " 3r « (0 '** ) + (e *(°) ' fa ( r )) 2 ] + 2r a(0,t*) r a 



(37) 



3r^(0, f) + 3 (e (0) • 2„ (**))<« - 4 (e (0) • a?„(t*)K 



c r 



r a -3 (e (0) ■ x a (t*)) 



C 2 T 2 + 2 C 3 T 3 



(38) 



Here, the sufBx a labels the masses and r a is the distance between the position of the photon along its unperturbed 
trajectory and the position of the mass m a at the time t* . Explicitly the distance r a is given by 



r a = r a (r,t*) 



c t + e - 2cre (0) • f„(t*) - 2£ • x (0 + a£(t*) 



1/2 



(39) 



It follows from the expression for r a that r a (0,t*) is the distance between the point of closest approach of the 
unperturbed light ray to the origin of the coordinate system and the position of the mass m a at the time t* . The 
quantity R a appearing in Eqs 1(55 )1 -1(55 )1 is given by 



R a =r 2 a (0,t*)-(e (oy x a (t*)) 2 



(40) 



To get the expression for the angle of light deflection we have to introduce 5Z?x)i( r ) an< ^ com P u te the 

limit t — ► oo. We have to compute the limit r — > oo, because we assume that the observer is located at infinity in an 
asymptotically flat space. 

After introducing the perturbation SI^^t) into Eq. (|18f) and computing the limit for r — > oo, we find 



a— 1 a 

2 

+ 4 E tt^) • - P X(**)] 

a— 1 

- 4 E f^c*) [* - p»)] - 4 E ^(g(o) • *un) 2 [f - pixun) 

a— 1 a— 1 



+ 4 E ir( e >) • «»(**)) 2 [e - pfrun] 

a— 1 

a— 1 " 

+ 4 E S^O) • Za(**))(e ( 0) • [f ■ t? (0 - *„(«*) • V a (t*)] [C - P*xi(fj\ 

a— 1 " 
2 

+ 4 E S • ^^)) 2 - Va(t*)(Xa(n ■ Va(t*)) + (*«(**) > 

a— 1 u 



(41) 
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V. THE POST-LINEAR LIGHT DEFLECTION IN THE POST-LINEAR GRAVITATIONAL FIELD OF 

TWO BOUNDED POINT-LIKE MASSES 

In this section we present our computations for light deflection in the post-linear gravitational field of two bounded 
point-like masses. In our computations we assume that both the light source and the observer are at infinity in 
an asymptotically flat space so that the effects of h^u and h$ near the light source and near the observer are 
negligible. We take into account only the terms of the order G 2 /c 4 . From equations Ijltjfl . (|17fl and (|18fl we see that 

a^ 2 -j is a function of the post- linear metric coefficients h^u and of the linear metric coefficients h^J . To facilitate the 
computations we separate the light deflection terms that are functions of the post-linear metric coefficients from the 
terms that are functions of the linear metric coefficients and the perturbations of the first order in G of the vector 
tangent to the unperturbed light ray. First we compute the terms of a 1 ^ that are functions of the post-linear metric 

coefficients, which we denote by aL r 



A. The post-linear light deflection terms that depend on the metric coefficients quadratic in G 

It follows from equations (|16fl and (|18|) that a part of the post-linear light deflection is given by: 



«(2)I 



-PI. 



drh. 



(2) I , 

00,g k-o + 



dr 



-M 2) - /i (2) 

2 mn : q qm,n 



i m in I 

W(o) k-o 



(42) 



Upon introducing the post-linear metric (|19fl and l|2(J|) into 142f) we obtain integrals whose integrands are functions of 
the distances r\, f2, S and their inverses. Through the distances ri, r-i and 5, the resulting integrals from Eq. I|42|) 
are functions of the positions of the masses x a {t). 

For the same reason as in the case of the linear light deflection we are here allowed to fix the values of the positions 
of the masses x a (t) to their values at the time t* before performing the integration. The resulting integrals are given 
explicitly in Appendix lAl 

To evaluate the integrals that cannot be represented by elementary functions we resort as usual to a series expansion 
of the integrands. To perform the series expansion we consider the integrands as functions of the distances r±, T2 
and S. Then, we expand these functions in a Taylor series about the origin of the coordinate system x\ — x*2 — to 
the second order. We need only to perform the Taylor expansion up to second order, since with an expansion to this 
order we obtain a result which is sufficiently accurate for the applications that we shall consider in this paper. 

The positions of the masses in the centre of mass frame defined in Subsection IIII Bl are given by equations i|24[l and 
(I25|l . Here, we do not need to take into account the lPN-corrections in the positions of the masses, because if we 
introduce these into Eq. I|42|l we will obtain terms of higher order than G 2 /c 4 . 

In Section I Villi we shall compute the post-linear light deflection terms resulting from the introduction of the 
lPN-corrections in the positions of the masses into the equation for the linear light deflection. 

Also, we do not need to consider here the correction terms arising from introducing the motions of the masses into 
Eq. H42f> . because these terms are of higher order than G 2 /c 4 . The correction to the post-linear light deflection arising 
from introducing the motion of the masses into the expression for the linear perturbation is denoted by a^ 2 j m and 
we shall compute it in Section IVTl 



B. The post-linear light deflection terms that are dependent on the metric coefficients linear in G 

The post- linear light deflection terms which are functions of the linear metric coefficients and the linear perturbations 
51(i)(t) we denote by aLw r It follows from equations lfH)|) and (|T%|) that the resulting expression for the post-linear 
light deflection o^ 2 ^ n is given by 



1 



-PI 



aril "oo,ml(^) 



dr 
dr 



h^{h^ n --h^)]l^ 0) \^ ) 



! (1) 

l 0m,q 



L 0q,m 



^)(r)| M 



12 



l(-0 



^O(0)<"(i)MI(-) 



^■Op,m'(0)^(0)^a)( T )l(^) 



(43) 



To compute cc^n' we introduce the expressions for the perturbations SI^Jt) given by Eq. I|t)3l) and the metric 
functions JUJ, (|23|l and 11' I'll into the expression for Qf( 2 )n- Here, we may use the same approximations as before, i.e. 
we can fix the values of the positions and velocities of the masses to their values at the time t* before performing 
the integrals. The resulting integrals are given in Appendix El As explained in the preceding section, with the help 
of a Taylor expansion of the integrands we can evaluate the integrals, which cannot be represented by elementary 
functions. 



VI. LIGHT DEFLECTION AND THE MOTION OF THE MASSES 



In this section we compute the correction terms to the linear and post-linear light deflection arising from the effect 
of the motion of the masses on light propagation. The correction terms to the linear and post-linear light deflection 
can be found by means of Taylor expansions of the linear perturbation (|31|l in which the coefficients depend on the 
sources' coordinates x % a and their successive derivatives with regard to t, namely 

dx a i, d2 < dv l a 

taken at the time t* . 



A. The linear light deflection and the motion of the masses 
The correction terms to the linear perturbation arising from the Taylor expansion of Eq. Ij^ljl are given by 
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(e (0 ) • v a (t*)) |cct - (e (0 ) • x a {t*)) 
30 _ ^ 6 
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— (r„ • w a (t*))cr 



'(0) 



+ (e (0) - t ; a (r)) 2 
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-rrcr + -r(*a ■ v a (t*))a z 
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(e (0) • ff„(t*))c<7 + • #„(**) - x a {t) ■ v a {t*) 



~^{r a - v a (t*))a 
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\a + -|(f a • v a {t*))a 2 



l(-) 



G 



a=l 
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(0) 



--|(iW«(t*)) s 



G 



a=l 
15 



G 



15 _ 



-7(r a ■ tf a (t*)) 4 - -r(r a ■ v a {tlY<{e) + -rvtm 



(r a -v a (t*)) , /3(r Q -i/ Q (i*)) 2 2 



2 r 3 



(0) 



(44) 



Considering that in the present work we compute the post-linear light deflection up to the order G /c , we must 
retain the linear light deflection terms up to the order G/c 4 . Notice that the linear terms of the order G/c 4 are of 
the same order as the post-linear terms of the order G 2 /c 4 since for a system of bounded point-like masses, the virial 
theorem applies (i.e. v\ ~ G/d) and, considering that the terms of the order G/c 4 are also terms in v 2 ., it is easy to 
see that these terms are of the same order as the post- linear terms of the order G 2 /c 4 . Linear light deflection terms 
of course we get from the perturbation ^muC 7 ") too. To obtain the perturbation S1 1 ^ u (t) we have to evaluate the 
integrals in the expression above. Taking into account that we are only interested in the angle of light deflection to 
the order G/c 4 , we need only retain the terms of the order G/c 2 and G/c 3 , since the expression for the light deflection 
angle (Eq. I|49|) ) contains a further factor 1/c. After performing the integration of the expression above and retaining 
only terms of the order G/c 2 and G/c 3 , we finally find 
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(e (0) • v a {t*))F a2 + £ • v a (t*) - x a {t*) ■ v a (t*) 
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(45) 



where the functions A a , F a i, F a2 , F a3 are given in Section llVl by Eqs The functions G a 2, G a 3 and G a 4 are 

given by 
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(46) 
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(47) 
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(48) 



After introducing the perturbation Sl^ u (r) into Eq. (|18f) and computing the limit for r — > oo, we obtain 
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(49) 



where the quantity R a is given by Eq. H4U|) . 



B. The post-linear light deflection and the motion of the masses 



The correction terms to the post-linear light deflection to the order G 2 /c 4 resulting from the Taylor expansion of 
Eq. ipIT|) are given by 
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(50) 



In the preceding equation the second integral diverges. Therefore we resort to a Taylor expansion of its integrand 
about the origin of the coordinate system x a — 0, up to second order. Then only the first term of the Taylor expansion 
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is a divergent integral and it is given by 
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Because in this case we do not need to take into account the lPN-corrections in the positions of the masses we can 
assume that the origin of the coordinate system which is located at the lPN-centre of mass coincides with the position 
of the Newtonian centre of mass. Taking into account the consequence of the Newtonian centre of mass theorem 



^ m a a a = 0, 



it is easy to see that the divergent integral vanishes. After performing the integration with respect to the parameter 
r we find: 



a ( 2 )m = A E ir^o) • Mt*Mm • s«(0) [e - p l A(t*)] 

a—1 
2 

- 4 ? E Ir^o) • f «(**)) 3 ( e >) • WW - p«<(n] 

a=l a 
2 

- 2 5 E 7T k"*' ~ • [ f - W'*)] 

a—1 
2 

" 4 § E if^o) • *«(**)) 2 [f- *.(**) - *»(**) • a a (t*)} [C - fyXt*)] 

a=l a 



G 



5E- a 2 - 



X 2 a (t*) , „(g(0)-^a(f)) 2 , (K(**)) 



e 2 



c 1 



(51) 



where the quantity R a is given by Eq. H4U|) . Note that in Eq. I|5U|I two of the three integrals were exactly integrated, 
so that the resulting expression given by Eq. I|51|l is a combination of exact terms with a term which is represented 
as an expansion in powers of (x a (i*)/£). In view of further applications and for the sake of uniformity we perform 
the Taylor expansion of the exact terms about the origin of the coordinate system x a = up to second order. After 
performing the Taylor expansion of the exact terms in Eq. ifBTll about the origin of the coordinate system x a = 0, up 
to second order, we find 
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+ 2 (52L^ + 12 (^))!} [? ._p. I ; (( . )] 

^t^d^-^^^^^pty^n (52, 

To replace the accelerations in the preceding equation by functions of the positions we shall use the Newtonian 
equations of motion. After replacing the accelerations in Eq. I|52l) and expressing the positions of the masses by their 
centre of mass frame coordinates without considering the lPN-corrections we find, 
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where the quantities r?i 2 , wi 2 and ri 2 are taken at the time t*. 
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VII. THE POST-LINEAR LIGHT DEFLECTION AND THE PERTURBED LIGHT RAY TRAJECTORY 

If we introduce into the equations for the linear perturbations 115|) and 14411 the expression for the perturbed light 
ray trajectory, we get additional post-linear light deflection terms. To compute these terms we have first to find 
the expression for the perturbation of the photon's trajectory that is linear in G. We obtain this perturbation by 
integrating the expression for the total linear perturbation given by Eq. (|63J) with respect to the parameter r. The 
expression for the total linear perturbation is obtained by summing up the expressions for the linear perturbations 
^(i)i ( r ) an d <^(i)n ( r )- Considering that in this paper we compute the post-linear light deflection to the order G 2 /c 4 
we do not need to retain in the expression resulting from the integration of Eq. I|tj3fl the terms of the order G/c 4 , 
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since these terms are related to the post-linear light deflection terms of higher order than G 2 /c 4 . After performing 
the integration of Eq. with regard to r and retaining only terms of the order 0(G/c 2 ) and 0(G/c 3 ) we obtain 
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where the functions B a , T a 2 and T a z are given by 
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As our integration constant we have chosen in Eq. 154|) 
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+ (l 

because with this integration constant we recover from our expression for the post-linear light deflection the correct 
expression for the post-linear light deflection in the event that the value of one of the masses is equal to zero (i.e. the 
Epstein-Shapiro post-linear light deflection). 

It follows from Eq. I|18|) that the expression for the linear light deflection for an observer located at infinity is given 

by 

4 )= ta {^^(r)}, 
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where the perturbations 5/^j(r) and SI%\ u (t) are given by Eqs iflS^ and 

After introducing the perturbation Sz^ into the equations for Sl^^ir) and SI^ u (t) we get a perturbed linear light 
deflection. Because the perturbation Sz^ is a small quantity compared to z(T) unpcrt ., we can resort to a Taylor 
expansion of the perturbed linear light deflection about Szm — in order to get the terms of the perturbed linear 
light deflection that are quadratic in G. We denote these terms by aj 2 , IV and they are given by 
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In the equation above we need not write out the expression G/c 3 {...} explicitly since it contributes only terms of 
order greater than G 2 /c 4 . After substituting the perturbation 154|) into the preceding equation we obtain the integrals 
for the post-linear light deflection a^ 2 ) Iv . Here, we take into account only the integrals of the order G 2 /c 4 . These 
integrals are given in an explicit form in Appendix IU1 



VIII. LIGHT DEFLECTION AND THE CENTRE OF MASS 



In this section we compute the corrections to the linear and the post-linear light deflection resulting from the 
introduction of the lPN-corrections to the positions of the masses in the equations for the linear perturbations 
(t) and <H/]\ n (T) given by Eqs l(T5)) and lfPE|) . It follows from equations and (|25T) that the lPN-corrections 
in the positions are given by 
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From expression l|58|) . it is easy to see that the corrections vanish when mi = m-i- The corrections also vanish for the 
case of circular orbits. After introducing the lPN-corrections into the equations for the linear perturbations SI^^t) 
and SI%^jj(t) we obtain the expression for the perturbed linear light deflection. Because the corrections 5x a are small 
quantities compared to x a , we can resort to a Taylor expansion of the perturbed linear light deflection about 8x a = 
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in order to find the correction terms for the linear and post-linear light deflection. We denote these terms by o^i^ 2 ) 
and they are given by 
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(59) 



For the same reason as in Eq. JSJJ, we do not need to write out the expression G/c 3 {...} explicitly in Eq. (|59|l . After 
substituting the lPN-corrections given by Eq. (|58|l into the equation above and taking into account only the terms of 
the orders G/c 4 and G 2 /c 4 we find, 
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Here, we have already replaced the photon trajectory by its unperturbed approximation z(r) unp ert. = t /(o) + £• 
After performing the integration we obtain, 
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where the quantity R a is given by Eq. i|40|) . Finally, considering further applications we express the positions of the 
masses by their centre of mass coordinates and expand the preceding expression about the origin of the coordinate 
system to the second order in (r^/c) to obtain, 
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As in the preceding section, the quantities ni 2 , V12 and ri 2 are taken at the time t* 



IX. THE TOTAL LINEAR PERTURBATION AND THE LINEAR LIGHT DEFLECTION IN THE 
GRAVITATIONAL FIELD OF TWO BOUNDED MASSES 

To obtain the total linear perturbation (r) we have to sum up the expressions for the linear perturbations 
^(i)i( T ) an< ^ ^(i)ii( r ) gi yen by Eqs l|32f> and (|45|l . The resulting expression is given by 
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In the expression above we need only retain the terms of the order G/c 2 and G/c 3 , since the terms of the order G/c 4 
are related to linear and post-linear light deflection terms of order higher than the terms which we compute in this 
paper. After introducing the perturbation SI^(t) into Eq. I|18(l and computing the limit for r — > oo, we find 
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where i? a is given by Eq. ljI0[l, 

Considering that the expression for the post-linear light deflection computed in this paper is given in terms of 
the centre-of- mass- frame coordinates and as an expansion in powers of (t" 12 /£), we must perform the expansion of 
the expression above about the origin of the coordinate system x a — and express the positions of the masses in 
terms of their centre of mass coordinates. After expressing the positions of the masses by their centre of mass frame 
coordinates without considering the lPN-corrections and expanding the expression given by Eq. I|64|l about the origin 
of the coordinate system to the third order in (rxz/Q we finally obtain, 
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where the quantities rti 2) U12 and ri2 are taken at the time t* . Notice that in Eq. I|65(l we chose the order of the 
expansion in an arbitrary manner in order to show the structure of the terms belonging to the linear light deflection. 
In concrete applications we have to choose the order of the expansion of the linear light deflection in accordance with 
the accuracy reached by the post-linear light deflection. 

To obtain the total expression for the angle of light deflection linear in G we have to add to the preceding equation 
the correction terms arising from the part of SA-^^y which is linear in G. The expression for cK-n^) * s gi ven by 
Eq. 



X. THE POST-LINEAR LIGHT DEFLECTION 



The final expression for the post-linear light deflection in the gravitational field of two bounded masses is obtained 
by summing up the parts of the light deflection, which are given in the preceding sections and in Appendices A to C. 
The final expression for the angle of light deflection quadratic in G to the first order in (7*12 /£) , in which the positions 
of the masses are expressed in the centre of mass frame coordinates, is given by 
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where o^ 2 ^ v is given by the part of (i.e. of Eq. (EH), which is quadratic in G. Here, it should be pointed out 

that some terms which belong to the post-linear light deflection given by Eq. I|66|l (e.g. the term 4 G mim ' 2 



arc 



related to terms in v\ 2 of the linear light deflection given by Eq. I|65|) through the virial theorem. In the next section 
we shall group these terms together before computing the limit r\ 2 — > in Eqs. Ij65(l and Ij66(l in order to remove the 
formal divergences. As in the case of the linear light deflection the order of the expansion in Eq. (|66|l was chosen 
in an arbitrary manner in order to show the structure of the terms belonging to the post-linear light deflection. In 
concrete applications the order of the expansion of the post-linear light deflection is to be chosen in accordance with 
the accuracy reached by the linear light deflection. 



XI. RESULTS 



The resulting expressions for the angle of light deflection linear and quadratic in G are given in Sections IIXI and IXl 
respectively. To study the important features of the derived results, let us consider special cases. 
The results take a particularly simple form for the following cases: 



A ) the value of one of the two masses is e> 



to zero 



If in Eqs l(6"5j) and (|rJfJ|l we put the value of one of the masses equal to zero (e.g. mi = M, m 2 = 0) we obtain 
expressions for the light deflection angle for a static point-like mass, 
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where in this case the ADM mass Madm is equal to the mass M. The deflection angle linear in G is the well known 
"Einstein angle" (see 0]). The light deflection angle given by Eq. I|68|l is the post-post- Newtonian light deflection for 
a point-like mass, which was obtained for the first time by Epstein and Shapiro and by other authors in 1980 (see [l2| 
and [i3j). 



B) the light deflection when ri2 — * 



In this subsection we are going to compute the limit of the expression for the linear and post-linear angle of light 
deflection (Eqs Ij65(l and l)66|l ) in the event that the distance r\ 2 between the components of the binary goes towards 
zero (i.e. r\ 2 — * 0). As we explained at the end of the preceding section, we have to group together the terms of 
Eqs. H62[> . H65|) and (|66() in an appropriate manner in order to remove the formal divergences. By inspection, it is 
clear to see that the remaining terms in Eqs. I|62|) . I|65|) and (|66|l are given by 
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After grouping together the terms that are related through the virial theorem in Eqs. (|69|l . (|70() and i|71|) we get 
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After simplifying Eqs i|72|) and i|73|l we finally obtain 

«(i) = - 4 4 + ~4~FT7 2 ( mi ~ m 2) 
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is the ADM mass of the system with [i = m\m 2 jM and M = mi + m 2 . In Eqs (|69() - (|74|l . the quantities ni 2 , v\ 2 and 
ri 2 are taken at the time t*. Note that the second term in Eq. (|74|l goes to zero when ri 2 — > since the expression 
in brackets remains finite. At the end, we recover the Einstein-angle (i.e. Eq. (|67|l ) with the ADM mass as given by 
Eq. l(75|l as well as the Epstein-Shapiro angle (i.e. Eq. JHH1))- 



C) the values of the two masses are equal and the light ray is originally orthogonal to the orbital plane of the binary 



In this case we choose M/2 for the value of the masses in Eqs (|65fl and (|66[1 and assume that the light ray is 
originally propagating orthogonal to the orbital plane of the two bounded masses (i.e. e(o) ■ n\ 2 — 0, e( ) • v\ 2 = 0). 
After introducing the ADM mass in the resulting expression for the angle of light deflection linear and quadratic in 
G and rearranging the terms, we finally find 
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where in this case the ADM mass is given by 

Madm = M 
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The expression for the linear light deflection given by Eq. I|77[) was expanded to the order (ri 2 /t;) 12 in order to reach 
the accuracy of the post-linear light deflection given in Eq. H78|) . Note that an expansion to the order (ri 2 /£) 12 in the 
linear part of the light deflection which is of the same accuracy as an expansion to the second order in (r\ 2 /£) in the 
post-linear part implies that (fi 2 /£) ~ (GA/adm/c^) 1 ^ 10 - ^ n -^qs (|77fl - l|79|) . the quantities ni 2 , Wi 2 and ri 2 are taken 
at the time t*. Note, that in this case the correction arising from the shift of the lPN-centre of mass with respect to 
the Newtonian centre of mass (see Eq. H62H ) vanishes. 



D) the values of the two masses are equal and the light ray is originally parallel to the orbital plane of the binary 

In this case we choose the value of the masses equal to M/2 in Eqs. (|65|l and (|66|l and assume that the light ray is 
originally propagating parallel to the orbital plane of the binary (i.e. • rii 2 = 0). After introducing the ADM mass 
as given by Eq. I|79|l and rearranging the terms, we finally find: 
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TABLE I: The parameters of PSR J0737-3039. 



Orbital period Pt (day) 


0.102251563(1) 


Eccentricity e 


0.087779(5) 


Total mass itia + rriB (M@) 


2.588(3) 


Mass ratio R = mA/mB 


1.069(6) 
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where, the quantities ni 2 , V12 and ri2 are taken at the time t*. 

In Eq. i jgDj l the components e|, -Pgnf 2 an d -Pg u i2 °f the linear light deflection were expanded to the order (ri 2 /£) 12 , 
(^12/0 7 an d (T12 /£) 7 respectively in order to reach the accuracy of the post-linear light deflection given in Eq. Q8ip. 
As in the preceding subsection here the correction arising from the shift of the lPN-centre of mass with respect to 
the Newtonian centre of mass (see Eq. (|62H ) vanishes. 

Finally, we apply the formulae for the angle of light deflection i|67|) - i|81|) to the double pulsar PSR J0737-3039. The 
parameters of the pulsar PSR J0737-3039 (e.g. see Q) are given in table [I] 

We compute the angle of light deflection for the cases when the distance between the two stars ri2 is maximal 
and minimal. In our computations we assume that the masses of the binary's components are equal, i.e. that the 
mass ratio R is equal to 1. For the impact parameter, we choose £ = hr\2- In order to compute the angle of light 
deflection we have first to calculate 7*12 and V\2- Note that we only use Newtonian relations, since the uncertainties 
in observational data, although small, are nonetheless greater than the corrections that post-Newtonian corrections 
would yield. 

To compute r\2 we have to calculate the semi-major axis of the elliptical orbit by means of the following equation 
(because of the low accuracy of the observational data, only Newtonian relations can be used) , 
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(82) 



where a denotes the semi-major axis and T the orbital period. The preceding equation follows from Kepler's third 
law (e.g. see |15|). The relationships between the the distances ri2 max , ?"i2min and the semi-major axis a are given 
by, 



l2max = a(l + e), 
l2mxn = a(l - e). 

We obtain the corresponding velocities to J , 12max and J , 12m i n from the equations given by 
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respectively. 

After introducing the parameters of PSR J0737-3039 into the preceding equations and into the formulae for the 
light deflection given in this section we find the results presented in table ITT1 

Table [n] shows the angles of light deflection computed by means of Eqs 1)78(1 and IjSljl corresponding to the Epstein- 
Shapiro angles of ai l m _ fl - ] = —1.55 • 10~ 6 arcsec and ot l ^ E _ s ^ — —2.20 • 10~ 6 



\E-S) 



arcsec. If we define the corrections to 



the Epstein-Shapiro angle that we calculated in this paper by 5a 
1. light ray originally orthogonal to the orbital plane 
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, we find: 



da 



(2)- 
«(2). 



= (-1.0 • 10~ 7 e\ + 5 ■ 10-* P % q n\ 2 ) arcsec for £ = 5 ri 2max , 
= (-1.5 • 10~ 7 el + 7- l0- & P l n\ 2 ) arcsec for £ = 5ri 2min . 
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TABLE II: The angles of light deflection linear and quadratic in G are given in arcsec. We denote by a^^gj and «(2)(e-S) 
the Einstein angle and the Epstein-Shapiro angle. For the light ray originally orthogonal to the orbital plane we assume that 
4 ■ ni2 = 1. For the light ray originally parallel to the orbital plane we assume that ejo) ■ ni2 = 1. 



ri2max = 9.56 


• 10 10 cm 


ri2min = 8.02 • 


10 10 cm 
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10 7 cm/s 
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■ 10 7 cm/s 
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^ a (2)|| = — ' 10 _7 e^ + 2 • 10" 8 P\n\ 2 arcsec for £ = 5ri 2ma x, 
cto| 2)|| = -1.6 • 10~ 7 e* + 3 • 10~ 8 P\n\ 2 arcsec for £ = 5r 12n 
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From the values of Sa l ^j_ and ^(2)11 we see that the corrections to the Epstein-Shapiro light deflection angle are 
greater for the case when the light ray is originally orthogonal to the orbital plane. Details related to the measurement 
of the corrections to the Epstein-Shapiro angle will be discussed in the next section. 



XII. DISCUSSION AND CONCLUSIONS 



The angle of light deflection in the post-linear gravitational field of two bounded point-like masses has been computed 
to the second order in G/c 2 . Both the light source and the observer were assumed to be located at infinity in an 
asymptotically flat space. The light deflection linear in G has been exactly computed. It was shown that the expression 
obtained for the linear light deflection is fully equivalent to the expression given by Kopeikin and Schafer in |5| in the 
event that the velocities of the masses are small with respect to the velocity of light and the retarded times in the 
expression of Kopeikin and Schafer are close to the time of closest approach of the unperturbed light ray to the origin 
of the coordinate system. To evaluate the integrals related to the light deflection quadratic in G, which could not 
be integrated by means of elementary functions we resorted to a series expansion of the integrands. For this reason 
the resulting expressions for the angle of light deflection quadratic in G are only valid for the case when the distance 
between the two masses ri 2 is smaller than the impact parameter £ (i.e. r^/S, < 1). The final result is given as a 
power series in rn/^. The expression for the angle of light deflection in terms of the ADM mass to the order G 2 /c 4 
including a power expansion to the second order in (ri2/£), in which (7-12 /£) ~ (GMadm/c 2 ^) 1 / 10 is being assumed, 
is given in an explicit form for a binary with equal masses in the event that the light ray is originally orthognal to the 
orbital plane of the binary. The expression for the angle of light deflection in terms of the ADM mass to the same 
order is also given for a binary with equal masses in the event that the light ray is originally parallel to the orbital 
plane of the binary. For a light ray originally propagating orthogonal to the orbital plane of a binary with equal 
masses the deflection angle takes a particularly simple form. In the case when one of the masses is equal to zero, we 
obtain the "Einstein angle" and the "Shapiro-Epstein light deflection angle", as we do when r\2 — » 0. Application 
of the derived formulae for the deflection angle to the double pulsar PSR J0737-3039 has shown that the corrections 
to the "Einstein angle" are of the order 10~ 2 arcsec for the case when ri 2 /£ = 0.2, see table [H] The corrections 
to the "Epstein-Shapiro light deflection angle" lie between 10 -7 and 10 -8 arcsec, see table |H] We conclude that 
the corrections to the "Epstein-Shapiro light deflection angle" are beyond the sensitivity of the current astronomical 
interferometers. Nevertheless, taking into account that the interferometer for the planned mission LATOR [3 

will 

be able to measure light deflection angles of the order 10~ 8 arcsec, we believe that the corrections to the "Epstein- 
Shapiro light deflection angle" computed in the present work might be measured by space-borne interferometers in 
the foreseeable future. 
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APPENDIX A: THE POST-LINEAR LIGHT DEFLECTION a\ 2)I 



In the integrals that are given in this appendix as well as in Appendices [B] an d El we already replaced the photon 
trajectory by its unperturbed approximation z(t) — tZ(q) + £, where Z(q) is given by Zm) = ce*(o). The distances r± 
and r 2 are given by 

1/2 

+ e - 2 c t e (0) • x a (**) - 2 i ■ x a (t* ) + x 2 a (t* ) 

with a = 1 and a = 2 for the distances r\ and r 2 . The distance S is defined by 

S = ri + r 2 + r i2 . 

Here, ri2 = — a?2 (^*) I is the distance between the two masses mi and m 2 at the time t* and the unit vector 

ni2 is given by 

r1i 2 = — [f !(«*)- 3» (**)]■ 
ri2 

The positions of the masses in the centre of mass frame without considering the lPN-corrections are given by 

xi = X 2 ri 2 (t*) 

and 

x 2 = -X x r X2 (f). 

The integrals resulting from the introduction of the post-linear metric coefficients 1)190 and (120)1 into the expression 
for ckLxj given by Eq. (|^2*)l are: 
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APPENDIX B: THE POST-LINEAR LIGHT DEFLECTION a\ m 

After introducing the expressions for the perturbations Sl 1 ,^ (r) given by Eq. (|63() and the metric coefficients l|21|) 
and l(2"2")l into the expression for a 1 ^^ given by Eq. (f^jl we find, 
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The functions A%, A 2 , B>\ and B 2 are given in Section llVl 
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APPENDIX C: THE POST-LINEAR LIGHT DEFLECTION a 



(2)rv 



In this appendix we give the integrals of the order G 2 / c resulting from the introduction of the perturbation SzUs (r) 
given by Eq. (|54|) into the expression for a^mi §i ven by Eq. (|57|l . 
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where the functions Si and B 2 are given in Section IVIII 
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APPENDIX D: THE LINEAR LIGHT DEFLECTION TERMS ARISING FROM THE TERMS OF h, 
AND hi 1 ! WHICH CONTAIN THE ACCELERATIONS OF THE MASSES 



(i) 



As we mentioned in Section IlIII the terms of the metric coefficients and hn~a which contain the accelerations of 



'00 "P9 

the masses were introduced into the metric quadratic in G after substituting the accelerations by explicit functionals 
of coordinate positions of the masses by means of the Newtonian equations of motion. To get the light deflection 
terms arising from these terms in a form suitable for the comparison of our computations with the linear light 
deflection computed by Kopeikin and Schafer [jj, we compute here the light deflection resulting from these terms 
before performing the substitution of the accelerations. The terms of h^J and hf^q which contains the accelerations 
are given by 
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From Eq. Ijl8|l it follows that the linear light deflection is given by 
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where SI^Jt) is given by Eq. Ijl5|) . After introducing the metric coefficients (|D II into the equation above we obtain 
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As in Eq. I|5l)fl the second integral in the preceding equation diverges. After performing the Taylor expansion of the 
second integrand about the origin of the coordinate system x a = up to the second order and taking into account 
the Newtonian centre of mass theorem we perform the integration of Eq. (|D3(1 . As result we find 
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APPENDIX E: THE ANGLE OF LIGHT DEFLECTION IN THE LINEAR GRAVITATIONAL FIELD OF 

TWO ARBITRARILY MOVING POINT-LIKE MASSES 



In this appendix we show that the expression for the light deflection angle in the linear gravitational field of two 
arbritrarily moving point-like masses computed by Kopeikin and Schafer [f| is fully equivalent to the linear light 
deflection computed in this paper in the event that the velocities of the masses are small with respect to the velocity 
of light and that the retarded times in the Kopeikin and Schafer expression are close to the time of closest approach 
t* . The light deflection angle given in in the case that the source is located at infinity and the observer at z(t) is 
given by 
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(El) 



where w Q (s Q ) is the velocity of the a-th mass, is the unit vector tangent to the unperturbed light ray, r a {r, s a ) is 
given by r a (T, s a ) = z(t) — x a (s a ) and r a (r, s a ) is the Euclidean norm of r a (r, s a ). Here, s a is the retarded time for 
the a-th mass defined by the light cone equation for s a . The light-cone equation for s a is given by 

s a + r a (T,s a )=T + t*. (E2) 

From Eq. I|E1|) it follows that for an observer located at infinity, the angle of light deflection is given by, 
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where the quantity R a {s a ) is given by 



Ra(Sa) = rl(0, S a ) - (e (0 ) • X a {s a )) 2 
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Here, it is worthwhile to note that the preceding expression for the light deflection angle is equivalent to the expression 
given by Eq. (139) in Q. 

With the help of the light-cone equation l|E2|) and the relationship for the time of closest approach from [{| , 



t* = t--{S {0 yz{t)) 
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~ t - -e (0 ) • r a (t, s a ) - -e (0 ) ■ 'Xa(sa): (E5) 

it was shown that 

Sa - t* — -5(o) • X a (s a ) - — — 

c Zr a 

- -e( ) • x a (s a ), (E6) 
c 

where d a is given by d a = \z a (s a ) - x a (s a )\. 

If the velocities of the masses are small with respect to the velocity of light and the retarded times do not differ 
significantly from the time of closest approach t* , we are allowed to use the Taylor expansion of the quantity 

x* a (s a ) ~ <(**) + <(t*)( So - O + \ai(t*))(s a - t*f (E7) 



After substituting into Eq. l|E6j) . we find 

s a -t*~ i(e (0) • x a (t*)) + i(e (0) • v a (t*))(s a - f) + ^-(e (0) • a a (t*))(s a - t*f . (E8) 



Now we solve Eq. (|E8(I iteratively with respect to (s a — <*) to obtain 

^ - t* ~ ie (0) ■ s?„(t*) + i(e (0) • £«(**)) (£((>) • *?«(«*)) + o(^). (E9) 



After performing the Taylor expansion of the expression for the light deflection angle given by Eq. I|E3(I and taking 
into account Eq. (|E9(I . we finally obtain 
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It is easy to see that the expression above is equal to the sum of Eqs l|5ip. il64[l and l|D4|) . Note, that to obtain the 
term in P^et^(i*) of Eq. l)E10jl . it is better to sum up the corresponding terms in Eqs 1)50(1 and (|D3|) before performing 
the integration in order to remove the formal divergences. After summing up and performing the integration of these 
terms we get: 
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where the first term on the last line vanishes as consequence of the Newtonian centre of mass theorem. 
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